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1. INTRODUCTION 

Electrical power blackout has always been a major problem not only to the electrical power 
engineers, but also to the consumers. Power system blackouts will give massive effect not only to the 
operation of the electrical power system, but also to the day to day business that can lead to much worse 
problems such as economic losses. It has been found out that most of the electrical power blackouts that 
occurred at most part of the world were caused by voltage instability [1]-[3]. Voltage instability will take 
place due to the failure of the electrical power system to maintain the bus voltages remain the same after the 
electrical power system is being subjected to a disturbance [4]-[5]. For every operation of electrical power 
system, it is very important that the system operates within the allowable voltage range and as far as possible 
from the voltage critical point. If the power system reaches the voltage critical point, the voltage system will 
enter unstable region and hence, lead to blackout. However, due to the rapid development that occurred all 
over the world especially in developing countries, has resulted in the increasing of electricity load demand 
per year [6]. The incremental value of load demand especially reactive power (Q) load demand will make the 
electrical power system to operate close to the voltage instability limit [7]. This is where the importance of 
analysing voltage instability take place. Voltage instability analysis is very important in order to make sure 
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that the load buses in the power system operate within the allowable voltage range. A load bus is considered 
as a weak bus if it operates near the voltage critical point. Modal analysis technique is very useful to 
determine the weak load buses [8]. In modal analysis technique, the load bus participation factor is identified. 
The load bus with the highest participation factor is the weakest bus in the system and vice-versa. However, 
the modal analysis technique involes very detail and intensive measures and calculations such as the 
determination of reduced Jacobian matrices, determination of eiganvalues etc. To overcome this, this 
research will apply Artificial Neural Network (ANN) in the modal analysis technique. ANN will be used to 
determine the value of the load buses participation factors. 


2. RESEARCH METHOD 
2.1. Modal Analysis Technique 

The modal analysis technique was introduced by Gao, Morisson and Kundur in 1992 [9]. This 
technique intends to determine the participation factors of the load buses. The participation factors can be 
determined from the values of eigenvalue and eigenvector of the reduced Jacobian matrix. The higher the 
value of participation factor, the closer the load buses towards experiencing voltage instability. The 
procedures of the modal analysis technique are explained as follows: 


2.1.1. Reduced Jacobian Matrix (Jr) 

The first step in the modal analysis technique the determination of the reduced Jacobian Matrix (Jr). 
The Jacobian matrix that is obtainable form the Newton Raphson power flow method represents the injected 
real power (P) and reactive power (Q). The relationship is shown in Equation (1) [2], [10], [11]. 


Tes J 
ie 7 ie i ie (1) 


where 
AP is the incremental change in bus real power 
AQ is the incremental change in bus reactive power 
Ao is the incremental change in bus voltage angle 
AV is the incremental change in bus voltage magnitude 
In order to identify the reduced Jacobian matrix (Jr), the value of AP is being substitute with 0. 
Equation (2) explains it. 


Iso = is | ie (2) 


From Equation (2), the following Equation (3) and Equation (4) can be obtained. 

Ad = —Jps* Jey AV (3) 

AQ = JosA6 + Joy AV (4) 
Equation (5) below is formed by substituting Equation (3) into Equation (4) 

AQ = AVUav — Jas Jes Jpv] or AQ = J AV (5) 
where J, = [ov — Jos Ips" Jpy]. Equation (6) is formed by reorganizing Equation (5). 

AV = J, * AQ (6) 


Equation (6) shows the relationship between the incremental changes of voltage and reactive power. 
This proofs that the incremental value of reactive power will affect the value of voltages. 


2.1.2. Determination of the Most Critical Mode 

The second step in modal analysis is the determination of the most critical mode. The eigenvalues 
and eigenvectors of J, can be used to determine the modes of the power network. The lowest value of 
eigenvalue of J, determines the most critical mode of the power system [9], [12], [13]. Equation (7) depicts 
their relationship [2], [9]-[11]. 
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Jp = €An (7) 


where 

é is the right eigenvector of J, 

A is the diagonal eigenvalue of J, 
n is the left eigenvector of J,. 


2.1.3. Determination of the Participation Factor 

The final and main procedure of the modal analysis technique is the determination of the 
participation factors of the load buses. The bus participation factor is an indicator that shows the tendency of 
a particular bus towards voltage instability. It should be calculated at the bus that has the most critical mode. 
The bus participation factor is calculated by using Equation (8) [9],[11],[13]. 


Phi = 8: (8) 


where 

Pk; is the participation factor of bus k to mode i 
&; § the it" column right eigenvector of J, 

n; is the i" row of left eigenvector of J, 

In addition, the bus participation factor obtained from Equation (8) is shown in matrix form. The 
row of the matrix indicates the number of the bus. The column of the matrix shows the mode of the power 
network. The larger the value of bus participation factor of a bus, the closer the bus towards voltage 
instability and vice versa. 


2.2. Application of Artificial Neural Network (ANN) 

In this research, ANN will be used to predict the value of participation factor that has been 
determined from Sub-section 2.1. In order to achieve that, the ANN model must be trained to learn the 
relationship between the input quantities and the output [14]. For this purpose, approximately 500 training 
data were generated. The multilayer perceptron with backpropagation (MLPBP) ANN model is the most 
widely used of ANN model [15],[16]. In this research the MLPBP ANN model will be optimized in order to 
predict the participation factors. The inputs of the optimized MLPBP ANN are the vslues of real and reactive 
power of loads (P and Q loads). 

Figure | depicts the configuration MLPBP ANN network that is being made up of input layer (P and 
Q load in per unit), hidden layer and the output layer (values of participation factors). The prediction output 
values will then be being compared with the participation factors determined in Sub-section 2.1 in order to 
find the errors. The error is important in order to improve the ANN’s predictions values. MATLAB Neural 
Network Toolbox is used to run the optimized ANN model. The optimized ANN model in MATLAB is 
shown in Figure 2. In this research, the numbers of neuron in the hidden layer is 10 which is the default value 
in toolbox. 


2.3. IEEE 14-bus System 

The IEEE 14-bus system has been chosen as the test power system this paper. The modal analysis 
technique will be applied upon this test power system. IEEE 14-bus test power system consists of one slack 
bus (Bus 1), four voltage controlled buses (Bus 2, Bus 3, Bus 6 and Bus 8) and nine load buses (Bus 4, Bus 5, 
Bus 7 and Bus 9 until Bus 14) [17]. Figure 3 [17]-[20] depicts the diagram of this test power system. 


INPUT HIDDEN ™ ™ OUTPUT 


| LAYER TARGET VALUES 
| (PARTICIPATION 


FACTORS) 
7 \__ PREDICTION It 
SL OuTPUuT _ — 

| -———=S_ + 

Y ERROR 


Figure 1. Configuration of the MLPBP ANN network 
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Figure 3. IEEE 14-Bus Test Power System 


. RESULTS AND ANALYSIS 
3.1. Modal analysis result 


In this section, the results obtained from performing the modal analysis technique on the test bus 


system are presented in Figure 4. 
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Figure 4. Participation Factors for The IEEE 14-Bus Power System 


As can be seen in Figure 4, Bus 14, Bus 10 and Bus 9 have high participation factor with Bus 14 is 


the highest which is 0.3190. This means that Bus 14 has the highest tendency towards voltage instability. One 
more important information that can be learned from Figure 4 is that Bus 5 and Bus 4 are the two most stable 
load buses in this power system. It is also noticeable from the results presented in Figure 4 that only the load 
buses that possess the participation factor. This is because this modal analysis technique focuses on the 
relationship between the incremental changes of voltage and reactive power as stated in Equation (6). Since 
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the voltages of the slack and voltage-controlled buses are fixed prior to the load flow analysis, no 
participation factor is available on these buses. 


3.2. Prediction of Participation Factors by ANN 

In this paper, ANN is used to determine the participation factors as shown in Figure 4. The 
participation factors from both methods were compared and analysed. Figure 5 depicts the participation 
factors obtained from ANN and the comparison with the data in Figure 4. 
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Figure 5. Participation Factors Obtained from ANN and the Comparison with The Data in Figure 4 


As can be seen in Figure 5, ANN has the ability to predict the participation factors close to the 
participation factor in Figure 4. The ‘error’ in Figure 5 is the difference (absolute value) between the 
participation factor in Figure 4 with the values predicted from ANN. The highest value of error is only 
0.0009 recorded at Bus 5. This shows that this configuration of ANN can be used to predict the participation 
factors and hence, help in improving the voltage instability modal analysis technique. 


4. CONCLUSION 

The research has shown the importance of analyzing voltage instability in electric power system. 
The modal analysis technique is very valuable in determining which load bus is prone towards experiencing 
voltage instability. The results have also shown that ANN has the ability to predict load buses participation 
factors that can be used for determining weak buses. Once the weak load buses is determined, the engineers 
and power system planners can use the information to plan and manage the increased of load demand 
properly in order to avoid voltage instability from taking place. Voltage instability must be avoided from 
taking place because the consequences of voltage instability are fatal. The power systems blackouts that are 
caused by voltage instability might spread to the entire country. Aside from power distribution problems, the 
power systems blackouts that are caused by voltage instability will cause much worse such as economic 
losses. It is suggested in future research to consider the affect of distributed generator such as solar panels or 
wind turbine towards the stability of the voltage. 
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